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Abstract. - We study the connection between the parameters of the fractional Fokker-Planck 
equation, which is associated with the overdamped Langevin equation driven by noise with heavy- 
tailed increments, and the transition probability density of the noise generating process. Explicit 
expressions for these parameters are derived both for finite and infinite variance of the rescaled 
transition probability density. 



Introduction. — Heavy-tailed distributions, i.e., pro- 
bability distributions with power tails and infinite second 
moments, are an important tool for studying a number of 
physical, biological, economical and other systems whose 
behavior is determined by rare but large events [1-4]. In 
many cases the continuous-time dynamics of these sys- 
tems can be effectively described by the (dimensionless) 
overdamped Langevin equation 



x{t)^f{x{t),t)+at), 



(1) 



, where x{t) [x{Q) = 0] is a state parameter of the sys- 
' tem, f{x,t) is a deterministic function, and £^{t) is a ran- 
, dom noise defined by the infinitesimal increments A7y(f) = 
Jj*^^ dt'S^it') (t — > 0) that are assumed to be indepen- 
I dent on non-overlapping intervals and distributed with a 
' heavy-tailed distribution. Since the tails of these distribu- 
tions can not be neglected, the classical stochastic theory, 
which is based on the ordinary central limit theorem, is 
I not applicable to eq. (1). Specifically, if the increments 
are distributed according to a Levy stable distribution [5] , 
i.e., ^(t) is a Levy stable noise, then the probability density 
P{x, t) that x{t) = X satisfies the fractional Fokker-Planck 
equation [6-11] which can by written as 



/_!^da;e-*'=^/i(x) and T-'^{hk} = h{x) ^ (l/27r) x 
J_^dke^^'^hk represents the Fourier transforms, and a, 
[3 and 7 are the parameters of the stable distribution. 

Because of the generalized central limit theorem [13], 
the Levy stable distributions constitute an important but 
a particular class of heavy-tailed distributions. In this 
Letter, we show that the fractional Fokker-Planck equa- 
tion (2) is valid also for all noises ^(t) whose increments 
have heavy-tailed distributions. Explicit expressions for 
the parameters of eq. (2) are derived in terms of the 
asymptotic characteristics of these distributions. 

Definitions and basic equations. Our starting 
point is the generalized Fokker-Planck equation [14] 



d_ 

dt 



P{x,t) 



d_ 

dx 



f{x,t)P{x,t) + T-^{Pk{t)\nSk}, (3) 



-^f{x,t)P{x,t) 



Here, the Riesz derivative, d°'/d\x\°', is defined as [12] 
d''h{x)/d\x\°' = -J^-H|fc|"^fc}, a pair T{h{x)} = = 



which corresponds to the Langevin equation (1) driven by 
an arbitrary noise. The term "arbitrary" means that the 
independent increments A?7(jr) = r]{jT -\- t) ~ JjUt) = 
Jj^^^ dt'^{t') (r ^ 0, j = 0,1,...) of the discrete-time 

noise generating process rjlriT) — X]j=o ^vU''') — 
1,2,.. .) are distributed according to an arbitrary proba- 
bility density function p{Ar], r). In other words, p{Ar], r) 
is the transition probability density of the process rjlnr). 
It is assumed that (i) p{Ar], t) is properly normalized, i.e., 
d{Arj) p{Ari , t) — 1, (ii) the first moment, if it ex- 
ists, equals zero, i.e., J^^d{Ari) p{Ari,T)Ari — 0, and 
(iii) limr^o p{Ari,T) — (5(Ary), where S{-) is the Dirac S 
function. The characteristic function Sk = [e~'^^^^^^) of 
77(1) — hnv^o X]j=o'~^ ^^(j''') ([!/■''] denotes the integer 
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part of 1/r) is connected with the characteristic function 

Pk{T) = {e~^'^^^^^^^) of Arj{jT), i.e., the Fourier transform 
of p{Ar], r), via the relation [14] 



InS'fc = lim - \pk{T) - 1]. 

T— »0 T 



(4) 



The generalized central limit theorem [13] implies that 
for a wide class of properly scaled probability densities 
^(Ar;, r) the characteristic function Sk corresponds to 
Levy stable distributions. These distributions are de- 
scribed by four parameters [5]: an index of stability a G 
(0, 2], a skewness parameter /3 G [—1, 1], a scale parameter 
7 € (0, oo), and a location parameter p G (— oo, oo), which, 
in accordance with the initial condition P{x, 0) = S{x), is 
assumed to be zero. Therefore, excluding from consider- 
ation the singular case when a = 1 and /3 ^ simulta- 
neously (in this case \ ln Sk\ = oo if k ^ 0, see below, 
and as a consequence the system reaches the final state 
immediately), one obtains 



Sk = exp 



-7|fcr (l + i/?sgn(fc)tan^)] , (5) 



and eq. (3) reduces to eq. (2) [14]. 

In order to find the stable parameters in eq. (2), the 
transition probability density p{Ari, r) must be specified. 
Next we consider a class of the functions p{Ar], r) defined 
as 



p(Ar?,T) 



1 



Arj 



(6) 



Here, a(T) is a positive scale function that vanishes at r = 
0, and the rescaled transition probability density, p{y), is 
an arbitrary probability density which satisfies the condi- 
tion lim(:^oP(y/e)/e = S{y) and has zero first moment (if 
it exists). According to this definition, }?fe(r) = Pka(T) and 
the normalization condition for p{y), which is equivalent 
to po = 1, yields InS'o = 0. lik^Q, then | ln5'fe| S [0, oo] 
and we can select three physically different situations de- 
pending on how quickly a{T) tends to zero as r ^ 0. First, 



if Pka( 



1 = o(r) (the scale function quickly vanishes). 



then InS'fc = and the noise is so weak that it does not 
effect the system at all. Second, if Pka{T) — 1 tends to zero 
slower than r (the scale function slowly vanishes), then 
llnS'fel = 00, i.e., the influence of the noise is so strong 
that the system relaxes instantaneously to the final state. 
Finally, the case we are primarily interested in corresponds 
to Pfea(r) ~ 1 = 0{t), i.e., < I InSfel < 00 and the noise 
acts on the system in a non-trivial way. 

Using eqs. (4), (6) and the representation In 5fe = R{k)+ 
il{k), we find the real. 



Rik) = - lim - 



dyp{y)[l-cos{ka{T)y)], (7) 



and imaginary, 



1 

I{k) = - lini - / dyp{y) sin(fca(r)j/). 



(8) 



parts of ln5fc. In the next two sections we will eval- 
uate R{k) and I{k) and express the stable parameters 
through a few main characteristics oi p{y) and a(r). It 
should be noted in this context that the results are quite 
different for the cases with finite and infinite variance 
= J^ao^yP^y^y^ ofthe probability density p(j/). Math- 
ematically, this difference arises from the fact that in the 
latter case the operations of taking the limit and integra- 
tion in eqs. (7) and (8) do not commute. 

Density functions with finite Vciriance. — If the 

variance is finite then, taking first the limit and then 
integrating, from eq. (7) we obtain 



R{k) = — 



lim 

r->0 



(9) 



Since, by assumption, the first moment of p(y) equals zero, 
for calculating J(fc) we temporarily assume that the third 
moment, TO3 = J^^dyp{y)y^, exists. This yields 



I{k) 



fc3 



lim . 



(10) 



As is seen from eqs. (9) and (10), the condition < 
I InS'fel < 00 holds only if a^(T) ~ (0 < g < 00). In this 
case I{k) = 0, T-^{Pkit)\nSk} = ia^q/2)d'^P{x,t)/dx'^, 
and the generalized Fokker-Planck equation (3) reduces to 
the ordinary Fokker-Planck equation [15] 

^^P{x, t) = t)P{x, t) + 7^ P{x, t), (11) 

which has the form of eq. (2) with 



a = 2, 7 = 



(12) 



If the third moment ms does not exist then I{k) can 
be evaluated by the method described in the next section, 
yielding the same result: I{k) = 0. Thus, if the incre- 
ments of the noise generating process have finite variance, 
the Langevin equation (1) is always associated with the 
ordinary Fokker-Planck equation (11). In this ease the 
noise £,{t) is white, i.e., it has a constant power spectral 
density at all frequencies, and 7 is the white noise in- 
tensity. In particular, for Gaussian white noise of inten- 
sity D characterized by the Gaussian probability density 
p{y) = (47r£')~^/^exp[-y^/(4£))] and the scale function 
a(r) = r^/^, we obtain cr^ = 2D, q = 1, and so 7 = £>. 

Density functions with infinite variance. — For 

noises with cr^ = 00 the power spectral density does not 
exist. To find R{k) and I{k) in this case, we consider a 
class of probability densities p{y) whose asymptotic be- 
havior is characterized by heavy tails, i.e., 



piy) 



{y ±00), 



(13) 
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with A = inin{A+,A_} G (0,2] and h± > 0. Let us first 
calculate R{k) at A € (0,2). In this case, representing 
eq. (7) in the form 



R{k) = - lim 



1 



dz 



air) 



+ P 



a{T)\k\ 



(1 — cos 2;), 



(14) 



replacing p{y) by the asymptotic formula (13), and using 
the integral relation 



dz- 



1 — cos z 



/o z^+^ 2r(l + A)sin(7rA/2) 

(A e (0,2), r(-) is the gamma function), we find 

^W- W 2r(l + A)sm(^A/2) r-S'o r ' 



(15) 



(16) 



where (5aa± is the Kronecker symbol. 

At A = 2 the integral in eq. (15) does not exist. There- 
fore, for calculating R{k) in this case, we present the inte- 
gral in eq. (14) as a sum of two integrals over the intervals 
(0,a(r)|fc|C) and (a(T)|A;|^, 00) with ^ = 0(1). This yields 
R{k) = Ri{k,(,) + R2ik,^), where 



and 



fc2 

I dy\p{y) +p{-y)]y^ \\m 
1 Jq r^o r 



2 a\r) 



(17) 



i?2(fc,0 = -fc2(/l+(52A+ + )lim^^ 

r— >0 T 



/>oo 

X / dz 

la(T)\k\£, 



1 — COS z 



(18) 



Since ..,.^dz{l — cos: 



(l/2)ln[l/a(r)] 



as r — *■ 0, the first term, Ri{k,^), can be neglected in 
comparison with the second, R2{k,^), yielding 

R{k) = 



,h+52A+ + h-62A_ a'^ir) , 1 
' + hm — — In 



r^O T 



a(r)- 



(19) 



In order to find explicit expressions for the imaginary 
part of InSk, we first rewrite eq. (8) in the form 



Ta{T)k Jo 



dz 



a{r)\k\ 



-P 



—z 



a{T)\k\ 



sm z. 



(20) 



Then, assuming that A G (0, 1), we substitute eq. (13) into 
eq. (20). Finally, taking into account that 



dz 



sm z 



z^+^ 2T{1 + A) cos(7rA/2) 



(21) 



if A G (0, 1), eq. (20) can be reduced to 



2r(l + A) cos(7rA/2) r^o r 



• (22) 



If A G (1,2] then the integral in eq. (21) diverges at 
the lower limit of integration, and the described approach 
becomes inapplicable. To generalize it to A G (1,2], we 
use the condition dyp{y)y = 0, which permits us to 
replace sin z by sin z — z in eq. (20). As a consequence, we 
arrive to the integral dz (sin 2; — z)z^^~^ that can be 
calculated by the same formula (21), i.e.. 



f 

Jo 



dz 



sm z — z 



^i+A 2T{1 + A) cos(7rA/2) ' 



(23) 



Thus, the representation (22) remains valid for A e (1,2] 
as well. 

At A = 1 both approaches developed for A G (0, 1) and 
A G (1,2] are not applicable (the integrals in eqs. (21) 
and (23) are divergent). Therefore, to calculate I{k) at 
A = 1, we use the method applied to find R{k) at A = 2. 
According to that we write /(/c) = /i(fc,^)+/2(fc,$), where 

h {k, = -k t dy \p{y) - p{-y)]y lim ^ (24) 



and 



h{k,C) = -k{h+6iA+-h_6iA-) lim 



T- Ja(r) 



dz 



smz 



(25) 

Then, using the asymptotic formula dz z ^ sin 2; ~ 

ln[l/o(r)] that occurs as r — »■ 0, one obtains 



I{k) = -k{h+6iA+ - h-6iA_ ) lim ^ In ^ 



r->0 T 



a(r) 



(26) 



Next, on the basis of the above derived results we can 
express the parameters of the fractional Fokker-Planck 
equation (2) through the asymptotic characteristics of the 
probability density p{y) (at y ±00) and the scale func- 
tion a(T) (at T 0). If A e (0,1) or A G (1,2) then, 
using eqs. (16) and (22) and the definition of the Riesz 
derivative according to which 



J^-'{\k\^Pk{t)} 



d\x\ 



-P{x,t), 



:F-'{^k\k\^-'P,m = -Q^Q^Pi^^t)^ (27) 

we obtain that the generalized Fokker-Planck equation (3) 
reduces to the fractional one (2) with 



a = A, P ■ 



h+SAA+ - /i-<5aa_ 

h+5 AA+ + h-S AA. ' 



7 = 



7r(/i+5 



AA 4 



h^5 



AA_ 



2r(l -I- A) sin(7rA/2) 



(28) 
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and q = limr^o o,^{t)/t. The condition < | InS'^l < oo 
assumes that < g < oo, and so the scale parameter must 
be proportional to r^^^, i.e., a(r) oc t^^^. We note also 
that p={h+- h_)/{h+ + h^) if A+ = A_ = A, /3 = 1 if 
A_ > A+ = A, and /3 = -1 if A+ > A„ = A. 

According to eqs. (16) and (26), at A = 1 and /i-|-5ia+ 
h-SiA- the condition < | In 5*^1 < oo holds only if a(T) oc 
T/ln(l/T). In this case R{k) = and eq. (3) takes the 
form 

5 . d 



dt 



P{x,t) 



dx 



[f{x,t) + fo]P{x,t), 



(29) 



where 



/o = g(/i+(5iA+ - h-SiA_ 



(30) 

and q = lim^^o a(r) ln[l/a(r)]/T. Equation (29) presents 

a very unexpected and truly remarkable result. It shows 
that, in contrast to common belief, a given type of noise 
does not generate a random dynamics of the system, but 
rather acts as a constant force /q. This force arises from 
the difference in the asymptotic behavior of the probability 
density p{y) at y ^ oo and y — > — oo. If /o = 0, i.e., 
= = h and A-|- = A_ = 1. then eq. (29) becomes 
invalid, and the next terms of the asymptotic expansion 
of p{y) as y ^ ±oo must be taken into account. For 
the special case of symmetric noise, when p{y) = p{—y), 
eqs. (8) and (16) yield I{k) = and R{k) = -\k\Trhq, 
respectively, and thus eq. (3) can also be written in the 
form of eq. (2) with 



a = 1, /3tan— =0, 7 = 7r/ig, 



(31) 



and q = lim^^o «(''')/''') i-e-, a{T) ~ qr. We note in this 
context that if the same scale function, a(r) ^ qr, is cho- 
sen for asymmetric noise with ^ /i_ then |/(fc)| — 00 
and, consequently, the system reaches the final state im- 
mediately. This result clarifies the nature of instabilities in 
the numerical simulation of x{t) which occur in the present 
case (see Ref. [16] and references therein). 

Finally, at A = 2 a non-trivial action of the noise (,{t) 
takes place only if the scale function satisfies the con- 
dition a^(T) ~ 2gT/ln(l/r) as r — *■ 0. In this case 
limT^oa^{T)\n[l/a{T)]/T = q, lim^-^Q a'^{T)/T = 0, and, 
as it follows from eqs. (19) and (22), R{k) — ~k'^{h-^^S2A+ + 
/i-52A_)g/2 and I{k) = 0. Therefore, at A = 2 eq. (3) 
takes the form of the ordinary Fokker-Planck equation 
(11), i.e., eq. (2) with 



2, 7 



h+S2A+ + h-S2A. 



(32) 



It should be emphasized that although eq. (11) is the same 
for 00 and A = 2, these two cases are quite differ- 

ent because at A = 2 the variance of p{y) is infinite. 
This difference results in different dependence of the scale 
parameter 7 on p{y) . Namely, while in the former case it 
is proportional to [see eq. (12)], i.e., an integral char- 
acteristic of p{y), in the latter case the scale parameter is 
determined by the tails oi p{y) [see eq. (32)]. 



Thus, we have shown that (i) each noise whose incre- 
ments have a heavy-tailed distribution acts on the sys- 
tem the same (in the sense of the probability density 
P{x,t)) as a certain Levy stable noise, and (ii) the action 
of each Levy stable noise can be reproduced by differ- 
ent noises characterized by different distributions of the 
noise increments. We have determined the parameters 
of the fractional Fokker-Planck equation (2) that corre- 
sponds to the overdamped Langevin equation (1). If the 
transition probability density of the noise generating pro- 
cess is heavy tailed then these parameters are expressed 
through the characteristics of the tails. Otherwise, the 
fractional Fokker-Planck equation reduces to the ordinary 
one. These theoretical results seems to be especially im- 
portant for the simplification and validation of numerical 
simulations of the Langevin systems driven by noises with 
heavy-tailed distributions of the increments. 
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